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We present a relativistic procedure for the chiral expansion of the two-pion exchange component 
of the NN potential, which emphasizes the role of intermediate ttTV subamplitudes. The relationship 
between power counting in ttN and A^'A'^ processes is discussed and results are expressed directly 
in terms of observable subthreshold coefficients. Interactions are determined by one- and two-loop 
diagrams, involving pious, nucleons, and other degrees of freedom, frozen into empirical subthreshold 
coefficients. The full evaluation of these diagrams produces amplitudes containing many different 
loop integrals. Their simplification by means of relations among these integrals leads to a set of 
intermediate results. Subsequent truncation to 0{q'^) yields the relativistic potential, which depends 
on six loop integrals, representing bubble, triangle, crossed box, and box diagrams. The bubble and 
triangle integrals are the same as in ttN scattering and we have shown that they also determine the 
chiral structures of box and crossed box integrals. Relativistic threshold effects make our results 
to be not equivalent with those of the heavy baryon approach. Performing a formal expansion 
of our results in inverse powers of the nucleon mass, even in regions where this expansion is not 
valid, we recover most of the standard heavy baryon results. The main differences are due to the 
Goldberger-Treiman discrepancy and terms of 0{q^), possibly associated with the iteration of the 
one-pion exchange potential. 



I. INTRODUCTION 

A considerable refinement in the description of nuclear interactions has occurred in the last decade, 
due to the systematic use of chiral symmetry. As the non-Abelian character of QCD prevents low- 
energy calculations, one works with effective theories that mimic, as much as possible, the basic 
theory. In the case of nuclear processes, where interactions are dominated by the quarks u and 
d, these theories are required to be Poincare invariant and to have approximate SU{2) x SU{2) 
symmetry. The latter is broken by the small quark masses, which give rise to the pion mass at the 
effective level. 

In the 1960s, it became well established that the one-pion exchange potential (OPEP) provides a 
good description of A^'A^ interactions at large distances. When one moves inward, the next class of 
contributions corresponds to exchanges of two uncorrelated pions [1] and, until recently, there was 
no consensus in the literature as how to treat this component of the force. An important feature 
of the two-pion exchange potential (TPEP) is that it is closely related to the pion-nucleon (ttN) 
amplitude, a point stressed more than thirty-five years ago by Cottingham and Vinh Mau [2] . This 
idea allowed one to overcome the early difficulties associated with perturbation theory [3] and led to 
the construction of the successful Paris potential [4], where the intermediate part of the interaction 
is obtained by means of dispersion relations. This has the advantages of minimizing the number of 
unnecessary hypotheses and yielding model independent results, but it does not help in clarifying 
the role of different dynamical processes, which are always treated in bulk. 

Field theory provides an alternative framework for the evaluation of the TPEP. In this case, one 
uses a Lagrangian, involving the degrees of freedom one considers to be relevant, and calculates 
amplitudes using Feynman diagrams, which are subsequently transformed into a potential. An 
important contribution along this line was given in the early 1970s by Partovi and Lomon, who 
considered box and crossed box diagrams, using a Lagrangian containing just pions and nucleons 
with pseudoscalar (PS) coupling [5]. A study of the same diagrams using a pseudovector (PV) 
coupling was performed later by Zuilhof and Tjon [6]. The development of this line of research led 
to the Bonn model for the A'^A'^ interaction, which included many important degrees of freedom and 
proved to be effective in reproducing empirical data [7]. On the phenomenological side, accurate 



potentials also exist, which can reproduce low-energy observables employing parametrized forms of 
the two-pion exchange component [8]. 

Nowadays, it is widely acknowledged that chiral symmetry provides the best conceptual framework 
for the construction of nuclear potentials. The importance of this symmetry was pointed out in the 
1970s by Brown and Durso [9] and by Chemtob, Durso, and Riska [10], who stressed that it constrains 
the form of the intermediate ttN amplitude present in the TPEP. 

In the early 1990s, the works by Weinberg restating the role of chiral symmetry in nuclear inter- 
actions [11] were followed by an effort by Ordofiez and van Kolck [12] and other authors [13,14] to 
construct the TPEP in that framework. The symmetry was then realized by means of non linear 
Lagrangians containing only pions and nucleons. This minimal chiral TPEP is consistent with the 
requirements of chiral symmetry and reproduces, at the nuclear level, the well known cancellations 
present in the intermediate nN amplitude [15]. On the other hand, a Lagrangian containing just 
pions and nucleons could not describe experimental tvN data [16] and the corresponding potential 
missed even the scalar-isoscalar medium range attraction [14]. 

One needed other degrees of freedom. The A contributions were shown to improve predictions 
by Ordonez, Ray, and van Kolck [17] and other authors [18]. Empirical information about the low- 
energy nN amplitude is normally summarized by means of subthreshold coefBcients [16,19], which 
can be used either directly in the construction of the TPEP or to determine unknown coupling 
constants (LECs) in chiral Lagrangians. The inclusion of this information allowed satisfactory 
descriptions of the asymptotic NN data to be produced, with no need of free parameters [20-23]. 

As far Eis techniques for implementing the symmetry are concerned, recent calculations of the 
TPEP were performed using both heavy baryon chiral perturbation theory (HBChPT) and covari- 
ant Lagrangians. In the former case [12,17,21-24], one uses non relativistic effective Lagrangians, 
which include unknown counterterms, and amplitudes are derived in which loop and counterterm 
contributions are organized in well defined powers of a typical low-energy scale. In this approach, 
relativistic corrections required by precision have to be added separately [25]. 

QCD is a theory without formal ambiguities and the same should happen with effective theories 
designed to be used at the hadron level. In the case of nuclear interactions, this allows one to expect 
that the chiral TPEP should be unique, except for the iteration of the OPEP, which depends on the 
dynamical equation employed. 

In the meson sector, chiral perturbation is indeed unique and predictions at a given order are 
unambiguous. However, the problem becomes much more difflcult for systems containing baryons. 
At present, the uniqueness problem is under scrutiny and two competing calculation procedures 
are available based on either heavy baryon (HBChPT) or relativistic (RBChPT) techniques. If 
both approaches are correct, they should produce fully equivalent predictions for a given process. 
Descriptions of single nucleon properties were found to be consistent, provided the nucleon mass 
is used as the dimensional regularization scale [26]. In the case of ttN scattering, comparison of 
predictions became possible only recently, through the works of Fettes, Meifiner, and Steininger 
[27] (HBChPT) and Becher and Leutwyler [28,29] (RBChPT). Differences were found, associated 
with the fact that some classes of diagrams cannot be fully represented by the heavy baryon series. 
Discussions of the pros and cons of these techniques may be found in refs. [30,31]. 

In the NN problem, all perturbative calculations produced so far were based on HBChPT 
[12,17,18,21-25]. On the other hand, an indication exists that A'^A'^ results are approach depen- 
dent, for the large distance properties of the central potential were shown to be dominated by 
diagrams that cannot be expanded in the HB series [32]. The main motivation of the present work 
is to extend the discussion of the uniqueness of chiral predictions to the B = 2 sector. We do this by 
calculating covariantly the TPEP to order 0{q'^) and comparing our results with the HB potential 
at the same order. 

Our presentation is organized as follows. In section II we give the formal relations between the 
relativistic TPEP and the intermediate ttN amplitude, whose chiral structure is analyzed in section 
HI. We discuss how power counting in ttA is transferred to the TPEP in section IV and how it is 
reflected into subthreshold coeflicients in section V. The problem of the triangle diagram, in which 
heavy baryon and relativistic descriptions disagree, is briefly reviewed in section VI. In section VII 
we discuss the dynamical content of the potential and the properties of important loop integrals 
used to express it. Our TPEP is suited to Lippmann-Schwinger dynamics and, in appendix C, we 
review the subtraction of the OPEP iteration, needed to avoid double counting. The full TPEP, 
which represents an extension of our earlier work [14,20], is derived in appendix D. This potential is 
transformed using relations among integrals given in appendix E and a new form is given in appendix 
F, which is simpler by the neglect of short range contributions. The truncation of these results gives 
rise to our C(g*) invariant amplitudes and potential components, displayed in sections VIII and IX. 



In section X we compare our TPEP with the standard heavy baryon version, using expansions for 

loop integrals derived in appendix G. Conclusions are presented in section XI, whereas appendixes 
A and B deal with kinematics and relativistic loop integrals. 



II. TPEP - FORMALISM 



The TPEP is obtained from the T matrix Ttp, which describes the on-shell process 
N{pi) N(p2) N(p'i) N(p'2) and contains two intermediate pions, as represented in fig.l. In 
order to derive the corresponding potential, one goes to the center of mass frame and subtracts the 
iterated OPEP, so as to avoid double counting. The NN interaction is thus closely associated with 
the off-shell ttN amplitude. 




FIG. 1. Two-pion exchange amplitude. 



The coupling of the two-pion system to a nucleon is described by T, the amplitude for the process 
'K°-{k)N{p) ■K^{k')N{p'). It has the isospin structure 

Tha ~ 5abT^ + iebacTc T 

and the evaluation of fig.l yields 
with 



2! J (27r)4 [fc2 - /i2] \J.,2 _ ^2] 



where /i is the pion mass and the factor 1/2! accounts for the exchange symmetry of the intermediate 
pions. The integration variable is Q = ik' + fc)/2 and we also define q = {k' — k), t = q^, and 
'^i = {p'i +Pi)-Q/2m. Our kinematical variables are fully displayed in appendix A. 
For on-shell nucleons, the sub amplitudes may be written as 



T± = u{p') [A^+ Q S±] u{v) 



and the functions and are determined dynamically. An alternative possibility is 



u{p) 



with = + vB^. This second form tends to be more convenient when one is interested in the 
chiral content of the amplitudes. The information needed about the pion-nucleon sub amplitudes 
A^, B^, and may be found in the comprehensive review by Hohler [16] and in the recent chiral 
analysis by Becher and Leutwyler [29]. 

The intermediate tvN subamplitudes A^, B^, and depend on the variables k^, k'^, v, and t. 
For physical processes one has k'"^ = fc^ = v ~> H and t < 0. On the other hand, the conditions of 
integration in eq.(2.2) are such that the pions are off-shell and the main contributions come from the 
region v Q. Physical amplitudes cannot be directly employed in the evaluation of the TPEP and 
must be continued analytically to the region below threshold, by means of either dispersion relations 
or field theory. In both cases one should preserve the analytic structure of the -kN amplitude, which 
plays an important role in the TPEP. 

The relativistic spin structure of the TPEP is obtained by using eq.(2.5) into eq.(2.3) and one 
has, for each isospin channel. 
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and 



/[...]= fJ^ I 



The Lorentz structure of the integrals I is reaUzed in terms of the external quantities q, z, W, 
and g'^'^, defined in appendix A. Terms proportional to q do not contribute and we write 
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These expressions and the spinor identities (A20) and (A22) yield 
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In order to display the ordinary spin content of this amplitude, we go to the center of mass frame 
and use identities (A32)-(A35), which allow one to rewrite Ttp, without approximations, in terms 
of the (2 X 2) identity matrix and the operators^ 

f, 2 (1) (2) r~, 2 /o (1) - (2) - (1) (2)n 

^LS = i ((t'^^+o-'^^)-^ X 2/4 ,Q,Q = o-'^^-qxz o-'^^-qxz . 

The two-component momentum space amplitude in the center of mass (CM) is derived by dividing 
T by the factor (4Bm), present in the relativistic normalization, and introducing back the isospin 
coefficients as in eq.(2.2). We then have the decomposition 
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^We use here the notation and results from Partovi and Lomon [5], eqs.(4.26)-(4.28). 



,± _ + Tern. _ .± I ^SS ± I ± ^LS ± S^Q ± 

with r"*" = 3 and t~ = 2. Finally, the momentum space potential, denoted by is obtained by 
subtracting the iterated OPEP from this expression, so as to avoid double counting. 



III. INTERMEDIATE -kN AMPLITUDE 



The theoretical soundness of the TPEP relies heavily on the description adopted for the interme- 
diate TriV amplitude. In this work we employ the relativistic chiral representation produced by the 
Bern group and collaborators [28,29,33], which incorporates the correct analytic structure. For the 
sake of completeness, in this section we summarize some of their results. 

At low and intermediate energies, the nN amplitude is given by the nucleon pole contribution, 
superimposed to a smooth background. Chiral symmetry is realized differently in these two sectors 
and it is useful to disentangle the pseudovector Born term [pv) from a remainder (R). We then 
write 

The pv contribution involves two observables, namely, the nucleon mass m and the -kN coupling 
constant g, as prescribed by the Waxd-Takahashi identity [34]. In chiral perturbation theory, de- 
pending on the order one is working with, the calculation of these quantities may involve different 
numbers of loops and several coupling constants^. Nevertheless, at the end, results must be orga- 
nized in such a way as to reproduce the physical values of both m and g in [35]. Following 
Hohler [16] and the Bern group, [33,29] in their treatments of the Born term, we use the constant 
g in these equations, instead of (gA/fn)- The motivation for this choice is that the wN coupling 
constant is indeed the observable determined by the residue of the nucleon pole. We write 
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where s = {p+k)^ = {p' + k')^ and u = (p — k')^ = {p' — k)^. The arrows after the equations indicate 

their chiral orders, estimated by using s — rn^ ~ WQ and u—rn^ ~ —W-Q, with W — P1+P2 ~ p'i+p'2- 
When the relative sign between the s and u poles is negative, these contributions add up and we 
have [l/(s — m^) — l/(w — m^)] 0{q~^). On the other hand, when the relative sign is positive, 
the leading contributions cancel out and we obtain [l/(s - m^) -|- l/{u - m^)] 0{q°). 

In ChPT, the structure of the amplitudes involves both tree and loop contributions. The 
former can be read directly from the basic Lagrangians and correspond to polynomials in and t, 
with coefficients given by the renormalized LECs. The calculation of the latter is more complex and 
results may be expressed in terms of Feynman integrals. In the description of nN processes below 
threshold, it is useful to approximate these contributions by polynomials, using 



where Xr stands for D^, B'^/v, Dj^jv^ or S^. The values of the coefficients be determined 

empirically, by using dispersion relations in order to extrapolate physical scattering data to the 
subthreshold region [16,19]. As such, they acquire the status of observables and become a rather 
important source of information about the values of the LECs. 



^For instance, up to 0{q'^) receives contributions from tree graphs of £^^^ - - - £^^^ and one-loop graphs from £' 
expressed in terms of its bare coupling constants. 



The isospin odd subthreshold coefiicieiits include leading order contributions, which yield the 
predictions made by Weinberg [36] and Tomozawa [37] (WT) for TriV scattering lengths, given by 



DwT = jj^^ 0{q) , (3.7) 

^WT-l^^ 0{q') . (3.8) 

Sometime ago, we developed a chiral description of the TPEP based on the empirical values 
of the subthreshold coefficients, which could reproduce asymptotic NN data [20]. As we discuss 

in the sequence, that description has to be improved when one goes beyond 0{q^). In nuclear 
interactions, the ranges of the various processes are associated with the variable t and must be 
accurately described. In particular, the pion cloud of the nucleon gives rise to scalar and vector 
form factors [33], which correspond, in configuration space, to structures that extend well beyond 1 
fm [32]. On the other hand, the representation of an amplitude by means of a power series, as in 
eq.(3.6), amounts to a zero- range expansion, for its Fourier transform yields only 5 functions and its 
derivatives. So, this kind of representation is suited for large distances only. At shorter distances, 
the extension of the objects begins to appear. 




FIG. 2. Long range contributions to the scalar and vector form factors. 



In the work of Becher and Leutwyler [29] we can check that the only sources of A'^A'^ medium range 
effects are their diagrams k and Z, reproduced in figure 2, which contain two pions propagating in 
the t channel. Here we consider explicitly their full contributions and our amplitudes and 



are written as 














= D+rit) + 






(2) + [dto'^^ + dt^uh + dt^t' 


Bt- 


= B^r{t) + 








Dr-- 


= D-,{t) + 




[doo 


V + diQV^ + do-^vt] , 


Bn = 


- Bra 


.(i) + 


[1/(2/,^) + 6,7,; 


(0) ^ 





(3) 



In these expressions, the labels (n) outside the brackets indicate the presence of leading terms 
of 0{q^), whereas the label mr denotes the contribution from the medium range diagrams of fig.2. 

This decomposition implies the redefinition of some subthreshold coefficients, indicated by a bar 
over the appropriate symbol. Their explicit forms will be displayed in the sequence. 
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FIG. 3. Dynamical structure of the ©(g*) ttA amplitude; the blobs represent terms coming directly from the effective 
Lagrangians. 

The dynamical content of the 0{q'^) T^n amplitude derived in [29] is shown in fig.3 and our 
approximation in fig. 4. In the latter, the first two diagrams correspond to the direct and crossed 
PV Born amplitudes, with physical masses and coupling constants. The third one represents the 
contact interaction associated with the Weinberg- Tomozawa vertex, whereas the next two describe 
the medium range effects associated with the scalar and vector form factors. Finally, the last diagram 
summarizes the terms within square brackets in eqs.(3.9)-(3.12). 



FIG. 4. Dynamical content of the approximate TriV amplitude. 



IV. POWER COUNTING 



One begins the expansion of the TPEP to a given chiral order by recasting the explicitly covariant 
Ttp into the two-component form of eq.(2.16). This procedure involves no approximations and one 
finds, in the CM frame, 
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The potential to order 0(g") is determined by 0{q"), {tss^tr'^ts} ^ 0(g ) and 



This means that one needs T^^ -> ©(g"), {T^"^^^,!^'^],*,!^^^*} -» ©(g"-^) and 



{T^^jT^^^^^} ^ 0{q"^). We now discuss how the chiral powers in those functions are related 
with those in the basic ttN amplitude. This relationship involves a subtlety, associated with the fact 
that and contain chiral cancellations. 



A generic subamplitude Ixy given by the product of the corresponding ttN contributions and 
we have 

= /[•••]{ [x^.] K] [x^] [y^] [x^] [y.i] [xi] [Y^] } ■ 

The loop integral and the two pion propagators, as given by eq.(2.11), do not interfere with the 
counting of powers, since /[•••] 0{q^). The loop iutcgration is symmetric under the operation 
Q — » —Q, which gives rise to the exchange s <-» m in the Born terms. In the case of [X^]''^^ [^p^]*'^', 
one is allowed to use 
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within the integrand. For the specific components this yields 
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These results show that, inside the integral, and Bp^ cannot be always counted as 0{q^) and 
0{q-^), respectively For the products [X^^f^ [^a and [X^f^ [Ypi]^'\ one uses 
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and has D±„ -> 0{q) and ^ 0{q-''). Assuming [y±](^' ^ 0{q''), one gets 
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Finally, in the case of [A^h]'*' [^i^j'^^ one just adds the corresponding powers. 
In this work we consider the expansion of the potential to 0{q*) and need — » 0{q'^), 
0{q^), and {I'-gj^^ ,1^%^} 0{q°). This means that, in the intermediate 
nN amplitude, we must consider to O(g^) and to 0{q). 



(4.8) 



V. SUBTHRESHOLD COEFFICIENTS 

The polynomial parts of the amplitudes to order 0{q^), as given by eqs.(3.7)-(3.10), are 
determined by the subthreshold coefficients of ref . [29] , which we reproduce below 
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where the parameters d and di are the usual renormaUzed coupUng constants of the chiral La- 
grangians of order 2 and 3, respectively [26]. The terms within square brackets labelled (mr) in 
some of these results are due to the medium range diagrams shown in fig.2 and must be neglected^, 

because we already include their contributions in D^r and -B^r- The terms bearing the (WT) label 



must also be excluded, for they were explicitly considered in eqs.(3.9)-(3.12). 
the redefinition mentioned at the end of section III. 



This corresponds to 
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TABLE I. experimental values for the subthreshold coefficients and medium range (mr) contributions in " units; experi- 
mental results are taken from ref. [16]. 
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^In rcf. [29], the contribution of the triangle diagram to djp includes both short and medium range terms and only the latter 
must be excluded. 



The values of the subthreshold coefficients are determined from itN scattering data and, in a 

chiral expansion to 0{jf'), they arc used to fix the otherwise undetermined parameters Ci and di. 
In our formulation of the TPEP, wc bypass the use of these unknown parameters, for the redefined 
subthreshold coefflcients are already the dynamical ingredients that determine the strength of the 
various interactions. This allows the potential to be expressed directly in terms of observable 
quantities. 

In table I wc show the experimental values of the subthreshold coefficients determined in ref. [16] 
and the sum of {WT) and (mr) contributions. The redefined values are obtained by just subtracting 
the latter from the former''. It is worth noting that the values of djj and 5qj are compatible with 
zero. 

When writing the results for the TPEP, it is very convenient to display explicitly the chiral 
scales of the various contributions. With this purpose in mind, we will employ the dimensionless 
subthreshold constants defined in table II. 



TABLE II. Dimensionless subthreshold coefflcients. 
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*We use gA = 1.25, = 93 MeV, /x = 139.57 MeV and m = 938.28 MeV. 



VI. RELATIVISTIC AND HEAVY BARYON FORMULATIONS 



In this section we review briefly the relativistic formulation of baryon ChPT and its relationship 
with the widely used heavy baryon techniques. Chiral perturbation theory is a systematic expansion 
of low-oncrgy amplitudes in powers of momenta and quark masses, goncrically denoted by q. The 
chiral Lagrangian consists of a string of terms, labeled by its power in q. To a given order, one 
builds the most general Lagrangian, consistent with Poincare invariance and other symmetries of 
QCD (parity, time reversal, and approximate chiral symmetry). A Lagrangian of order n produces 
tree graphs of the same order, while loop graphs are expected to contribute at higher orders, following 
a power counting scheme. This is indeed what happens in the mesonic sector, where loop graphs 
are two orders higher than tree graphs, if one uses dimensional regularization. 

In relativistic baryon ChPT, dimensional regularization no longer leads to a well defined power 
counting [33], loops start at the same order as tree graphs and the connection between loop and 
momentum expansion is lost. A similar phenomenon is observed in the mesonic sector if one uses 
another regularization scheme, such as Pauli-Villars. 

In HBChPT, this problem is overcome by means of the expansion of the original Lagrangian 
around the infinite nucleon mass limit [38]. One integrates out the heavy degrees of freedom of 
the nucleon field, eliminates its mass m from the propagator, and expands the resulting vertices in 
powers of 1/m. This formulation gives rise to a power counting scheme, but Lorentz invariance is no 
longer explicit. It can still be recovered, but only after a rcsummation of all terms in this expansion. 

The HB approach also has a more serious problem, pointed out recently by Becher and 
Leutwyler [28], namely, that it fails to converge in part of the low-energy region. In order to 
avoid this, they proposed a new regularization scheme, the so called Infrared Regularization, which 
is manifestly Lorentz invariant and gives rise to a power counting. The method is based on a previ- 
ous work by Ellis and Tang [31], where a loop integral H was separated into "soft", infrared (7) and 
"hard" , regular {R) pieces. The former satisfies a power counting rule and has the same analytic 
structure as H in the low-energy domain. The latter may contain singularities only at high energies 
— in the low-energy region, it is well behaved and can be expanded in a Taylor series, resulting in 
polynomials of the generic momentum q. Therefore the hard pieces, which are the power counting 
violating terms, can be absorbed in the appropriate coupling constants of the Lagrangian and one 
considers only 7, the infrared-regularized part of 77*. 

Ellis and Tang have shown that the chiral expansion of the infrared regularized one-loop integral 7, 
with the ratio q/ ji fixed, reproduces formally the corresponding terms in the HBChPT approach [31], 
even in the cases where such an expansion is not permitted. This allows one to assess the domain 
of validity of the HB series. 

For the sake of completeness, in the sequence, we reproduce some of the results derived by Becher 
and Leutwyler. They have analyzed in detail the triangle graph of fig. 4, which contributes to the 
nucleon scalar form factor, and shown that the HBChPT formulation is not suited for the low 
energy-region, near t = . Its exact spectral representation is given by [33] 
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*This problem has been recently reviewed by Meifiner in sects. 3.4-3.7 of ref. [30]. 



seems to be of order q~^, and the HB chiral expansion of (6.2) would yield tan~^ x = tt/2 — 1/x + 

l/Si^ + • • •. However, this representation of tan~^ x is valid only in the domain |a::| > 1. For |a;| < 1, 
one should use tan^^ x — x — x^/3 + • • •, but this corresponds to an expansion in inverse powers of 
q. From (6.3) wo sec that the HB expansion of (6.1) breaks down when t' approaches 
Becher and Leutwyler have shown that it is possible to write accurately 
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By keeping only the first bracket in the integrand, one recovers the heavy baryon result. However, 
the region t ~ 4/i'^ is dominated by the lower end of integration in f', where the second term becomes 
important. The HB approximation is not valid there. The integration can be performed analytically 
and Becher and Leutwyler found 
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with T = t/ jjb^ . This result is interesting because it shows clearly that, for values of t far from Ajjb^ , 
the contributions of the two brackets decouple and can be expanded in powers of q. The second 

term is then O(q^). On the other hand, when t ^ 4fi^, both contributions merge, the full result for 
7(t) is the outcome of large cancellations between them, and an expansion in q does not apply. In 
fig.5, we display the behavior of the various terms in eq.(6.5) in the range 3.5/i^ < i < where 
the second bracket is important. In this figure we also show the effect of making 



m^/T 



2 m 



- In 1 



+ ln 



m2(4 ■ 



4m2 



This rough approximation is not mathematically precise, but it allows one to guess the order of 
magnitude of the threshold contribution. 



(6.5) 
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FIG. 5. Behavior of the function 7(t) as given by eq.(6.5) (full line) and partial contributions: HB (dashed line), th (dotted 
line) and eq.(6.6) (dot-dashed line). 

The discussion of the behavior of the triangle diagram in the neighborhood of t = is relevant to 
the A^A^ potential because, in configuration space, this region describes its long distance properties, 
as observed numerically in our previous works [20,32]. To see this, let us take the representation of 
(6.1) in configuration space: 

^ ^ r f -iQ r Im7(f') 1 r ^"'^ t u'^ la 7^ 



The exponential in the integrand shows clearly that, for large values of r, results are dominated 

by the lower end of the integration. Thus, if we want to have a good description of r(r) at large 
distances, we need a decent representation for lm'y{t') near t' = which is not provided by 
HBChPT. 



VII. DYNAMICS 

The chiral two-pion exchange potential is determined by the processes depicted in fig.6, derived 
from the basic tvN subamplitudc and organized into three different families. The first one corre- 
sponds to the minimal realization of chiral symmetry [14], includes the subtraction of the iterated 
OPEP, and involves only pion-nucleon interactions with a single loop, associated with the constants 
m, g, and The same constants also determine the two-loop processes of the second family. The 
last family includes chiral corrections associated with subthreshold coefficients and LECs, represent- 
ing either higher order processes or other degrees of freedom. 
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(a) (b) (c) (d) (e) (f) 



(g) (h) (i) Q) 



'ii — i — V ^# 

' ' ' ' ^p^ 

(k) (1) (m) (n) (0) 

FIG. 6. Dynamical structure of the TPEP. The first two diagrams correspond to the products of Born nN amplitudes, 
the third one represents the iteration of the OPEP, whereas the next three involve contact interactions associated with the 
Weinberg- Tomozawa vertex. The diagrams on the second line describe medium range effects associated with scalar and vector 
form factors. The remaining interactions are triangles and bubbles involving subthreshold coefficients. 

The first two diagrams of fig.6, known, respectively, as crossed box and box, come from the products 
of the nN PV Born amplitudes, given by eqs.(3.2)-(3.5) and involve the propagations of two pions 
and two nucleons. The third one represents the iteration of the OPEP and gives rise to an amplitude 
denoted by Tu, derived after the work of Partovi and Lomon [5] and discussed in detail in appendix 
C. The remaining interactions correspond to triangle and bubble diagrams, which contain a single 
or no nucleon propagators, besides those of two pions. 

The construction of the TPEP begins with the determination of the relativistic profile functions, 
eqs.(2.7)-(2.10), using the nN subamplitudes and discussed in section III. Results are then 
expressed in terms of the one- loop Feynman integrals presented in appendixes B and C, which may 
involve two, three, or four propagators. The evaluation and manipulation of these integrals represent 
an important aspect of the present work and it is worth discussing the notation employed. 

Momentum space integrals are generally denoted by 11 and labeled in such a way as to recall 
their dynamical origins. We use lower labels, corresponding to nucleons 1 and 2, with the following 
meanings: c contact interaction; s s-channcl nucleon propagation; and u u-channel nucleon 
propagation. This means that functions carrying the subscripts (cc), (sc), (ss), and (us) correspond, 
respectively, to bubble, triangle, crossed box, and box diagrams. The last class of integrals includes 



the OPEP cut, which needs to be subtracted. This subtraction is implemented by replacing the 

(us) integrals by regular ones, represented by the subscript (reg) and given in appendix C. Upper 
labels, on the other hand, indicate the rank of the integral in the external kinematical variables q, z 
and W. For instance, the rank 2 crossed box integral is written as 

/ "'<3 f Q^Q'^ 1 1 2m/i 2mM 
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All integrals are dimensionless and include suitable powers of pion and nucleon masses, so as to 
make them relatively stable upon wide variations of the latter. We have studied these integrals 
numerically and, typically, they change by 30% when one moves the nucleon mass from its empirical 
value to infinity. The fact that the integrals are is rather useful in discussing chiral scales 

and heavy baryon limits. At present the infrared regularization techniques axe still being developed 
for the case of two nucleon system [39] and we have used dimensional regularization whenever 
appropriate. As a consequence, our results are accurate only for distances larger than a typical 
radius. Our numerical studies in configuration space indicate that this radius is of about 1 fm. 

The covariantly expanded TPEP, to be given in section X, is expressed in terms of the functions 
n^c""^ , n^c""^ , nis""^ , n^eg"^ , and niVf^ ■ in order to simplify the notation, in the main text we call 
them Hi, lit, Hx, Hi,, and lit, respectively. 

The function He represents the bubble diagram and is given by 
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This integral can be performed analitically® and its regular part may be written as 



n, = -2 ^] In f yr~^7v + v-w) ■ 

The function lit, associated with the triangle diagram, is expressed by 
f d*Q 1 2mM ^ i 



and related to the function '){t) discussed in the preceding section by lit = — 2m/u(47r)^7(i). The 
heavy-baryon representation of this function is 



nt-nr=n„ + -^Lnf^ 

2m 



with'^ 



= , tan-1 , 

iC = (i-V2m') m , 

and Yi' ^ ^l (dU/d^i). 

The functions IIx , 11;,, and llf, are associated with crossed box and box diagrams and their complete 
expressions are given in appendix B. Their heavy baryon expansions are derived in appendix G and 
read 



The function n« is related to the L{q) used in ref. [21] by Ue = -2L{q) and to the J{t) of [29] by Ue = {4tv)'^ J 
The function Ha is related to the A{q) of ref. [21] by Ha = — 47r/uA(g). 
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In the heavy baryon expansion of the potential, the following results are useful 

n^ = 2 + n^/(i-vV), 

n'l = 2/(1-V4m') + [2/(l-t/4M') - l/{l-t/4.u?f] , 

nl = n„-7ri/(i-i/4M^) . 

For the reasons discussed in the preceding section, all these heavy baryon representations are 
inaccurate around t ~ Ajx^ . 

VIII. COVARIANT AMPLITUDES 
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The direct reading of the Feynman diagrams of fig.6 gives rise to our full results for the relativistic 

profile functions, displayed in appendix D. These are the functions that the chiral expansion must 
converge to and hence they allow one to assess the series directly. On the other hand, they do not 
exhibit explicitly the chiral scales of the various components of the potential, since their net values 
are the outcome of several cancellations. 

In order to display these scales, in appendix E we derive several relations among integrals, which 
are used to transform the full results of appendix D into the forms listed in appendix F. The relations 
given in appendix E are, in principle, exact, provided one keeps short range integrals that contain a 
single or no pion propagators. However, for the sake of simplicity, we neglect those contributions*. 
The importance of this approximation was checked by comparing numerically the Fourier transforms 
of the various amplitudes of appendixes D and F. In all cases, agreement is much better than 1% 
for distances larger than 1 fm, except for T^jj, whore the difference is 4% at 1.5 fm and falls below 
1% beyond 2.5 fm. This has very little influence over the full potential. 

With the purpose of allowing comparison with results produced in the HB tradition, we write 
our final expressions for the potential in terms of the axial constant qa, which is related to the 
nN coupling constant by g = (1 + Agt) gAm/f-K- Here Agt is the Goldberger-Treiman (GT) 
discrepancy®, proportional to ji^ . In applications, on the other hand, we recommend the direct use 
of the -kN coupling constant g, by making gA = gfTr/m and neglecting Agt in our results. 

The appropriate truncation of the expressions of appendix F, at the orders in q prescribed at the 
end of section IV, leads to the following results for the profile functions: 
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It would be very easy to keep those terms, but this would produce longer equations. 
®The GT discrepancy may be written [29] as Agt = -2di_%tJ^ j g + (^(q*)- 
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The results for the basic subamphtudes presented in this section are closely related to the under- 
lying nN dynamics and, in many cases, this relationship can be directly perceived in the final forms 
of our expressions. For instance, reorganizing the contributions proportional to lit in eq.(8.10), one 
has 
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The terms within the parentheses represent the contributions from fig.4, which read: (a) Born 

terms, proportional to g\; (b) Weinberg- Tomozawa term; (c) two- loop medium range interactions; 
(d) other degrees of freedom plus two-loop short range interactions. The organization of the last 
three terms may be better understood by noting that, around the point t = 0, the following expansion 
holds: {l—t/4fj?) lit — » — 7r-|-t 7r/6/i^, and the content of the parentheses of eq.(8.12) may be written 
as 
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This shows that the structure of eq.(3.12) is recovered, except for the medium range contribution, 
which is divided by a factor 2, chaxacteristic of the topology of Feynman diagrams. 



IX. TPEP 

Our final result for the relativistic 0{q'^) two-pion exchange potential is obtained by feeding the 
truncated covariant profile functions of the preceding section into eqs.(4.1)-(4.5). It is ready to be 
used as input in other calculations and is expressed in terms of five basic functions (section VII) and 
empirical subthreshold coefficients (section V). If one wishes, the latter may be traded by LECs, 
using the results of section V. The various components are listed below. 
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This potential is the main result of this work. If one keeps only terms up to order 0{(f), it 
coincides numerically with that derived earlier by us [20]. As far as 0(5*) terms are concerned, 
the only difference is due to the explicit treatment of medium range contributions. In our previous 

study we have shown that diagrams (k)-(o) of fig. 6 strongly dominate the potential. In the above 
expressions, those terms are represented by products of g\ by subthreshold coefficients. About 70% 
of the isoscalar potential comes from the term proportional to (5q^ + ij jj?^, which is related 
to the scalar form factor of the nucleon [32] , given by 



The leading contribution to then reads 
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As the scalar form factor represents the probing of the part of the nucleon mass associated with 

its pion cloud, the leading term of the A'^A'^ potential corresponds to a picture in which one of the 
nucleons, acting as a scalar source, disturbs the pion cloud of the other. A rather puzzling aspect 
of this problem is that the largest term in a ©(g^) potential is of 0{(i'). 



X. COMPARISON WITH HEAVY BARYON CALCULATIONS 

The relativistic potential of the preceding section involves five basic functions, representing loop 
integrals, and subthreshold coefficients. The latter can be rooxprossod in terms of LECs and explicit 
powers of /n/m, using the results of ref. [29], summarized in section V. The loop functions were 
derived by means of covariant techniques and one uses the results of section VII and appendix B. As 
discussed by Ellis and Tang [31] and in our section VI, if one forces an expansion of the relativistic 
functions in powers of ^/m, even in the regions where this expansion is not valid, one recovers 
formally the results of HBChPT. This procedure amounts to replacing the relativistic functions, 
which cover the neighborhood of the point t = 4/j^, by the heavy baryon series, which is not valid 
there. 



Performing such a replacement in the 0{q*) results of the preceding section, we find (inequivalent) 
expressions that coincide largely with those produced by moans of heavy baryon techniques. In order 
to allow comparison with ffBChPT calculations, in this section wc display the full /i/m expansion 
of our potential, without including terms duo to the common factor m/E. 

We reproduce below the results of refs. [21,24,25], which include relativistic corrections and were 
elaborated further by Entem and Machleidt [40] . The few terms that are only present in our potential 
are indicated by [•••]*: 
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15367r2/4m2 
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10247r2/6 



5i(32M' + 14q') + (4M'+q') 
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+ 
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XI. SUMMARY AND CONCLUSIONS 



We have presented a 0{q^) relativistic chiral expansion of the two-pion exchange component of 
the NN potential, based on that derived by Becher and Leutwyler [28,29] for elastic -wN scattering. 
The dynamical content of the potential is given by three families of diagrams, corresponding to 
the minimal realization of chiral symmetry, two-loop interactions in the t channel, and processes 
involving nN subthreshold coefflcients, which represent frozen degrees of freedom. 

The calculation begins with the full evaluation of these diagrams. Results are then projected into 
a relativistic spin basis and expressed in terms of many different loop integrals (appendix D). At 
this stage, the chiral structure of the problem is not yet evident. However, chiral scales emerge when 
these first amplitudes are simplified by means of relations among loop integrals. This gives rise to our 
intermediate results (appendix F) , which involve no truncations and preserve the numerical content 
of the various subamplitudes for distances larger than 1 fm. The truncation of these intermediate 
results to C(g^) yields directly the relativistic potential (section IX), which is ready to be used in 
momentum space calculations of NN obscrvables. 

Our treatment of the iVA'' interaction emphasizes the role of the intermediate nN subamplitudes 
and, in this sense, it is akin to that used in the Paris potential. We discuss how power countings in itN 
and NN processes axe related (section IV) and results are expressed directly in terms of observable 
subthreshold coefficients. The LECs d and di are implicitly kept within these coefflcients, grouped 
together with two-loop short range contributions. 

If the potential presented here were truncated at order 0(q^), one would recover numerically the 
results derived by us sometime ago [20] . However, processes involving two loops in the t channel do 
show up at ©(g*) and results begin to depart at this order. 

The dependence of the potential on the external variables is incorporated into five loop integrals, 
associated with bubble, triangle, crossed box, and box diagrams. The triangle integral is the same 
entering the scalar form factor of the nucleon and can be represented accurately by means of elemen- 
tary functions (section VII) and has the correct analytic behavior at the important point t = 
We have shown that this kind of representation can also be used to disclose the chiral structures 
of box and crossed box integrals (appendix G). The effects associated with the correct analytic 
structure of relativistic integrals are important because they dominate the long distance behavior 
of the potential. 

The expansion of the functions entering the relativistic potential in powers of /x/m is not mathe- 
matically defined around t — Nevertheless, in order to compare our results with those produced 
by means of HBChPT, we have assumed that such an expansion could be made for all low-energy 
values of t. This expansion then reproduces most of the standard HBChPT results. We find, how- 
ever, two systematic differences, apart from some minor scattered ones. The first one is due to the 
Goldberger-Treiman discrepancy. The other one concerns terms of 0{q^), whose origin is less cer- 
tain. However, the fact that they occur at the same order as the iteration of the OPEP suggests that 
there may be an important dependence on the procedure adopted for subtracting this contribution. 
This aspect of the problem is rather relevant in numerical applications of the potential and deserves 
being clarified. 

The numerical implications of the various approximations required to derive the 0{q'^) potential 
in configuration space will be presented in a forthcoming paper. 

ACKNOWLEDGEMENTS 

We thank C. A. da Rocha for supplying his numerical profile functions for the NN potential and 

J. L. Goity for useful discussions. R. H. also acknowledges helpful communications with T. Becher 
and M. Mojzis, the kind hospitality of the Theory Group of Thomas Jefi^crson National Accelerator 
Facility, and the financial support by FAPESP (Fundagao de Amparo a Pesquisa do Estado de 
Sao Paulo). This work was partially supported by DOE Contract No. DE-AG05-84ER40150 under 
which SURA operates the Thomas Jefferson National Accelerator Facility. 

APPENDIX A: KINEMATICS 

The initial and final nucleon momenta are denoted by p and p' , whereas k and k' are the momenta 
of the exchanged pions, as in fig.l. We define the variables 



W = pi + P2 = Pi + p'2 , 

Z = [{Pl + P'l) - {P2 + p'2)]/2 , 

q = k' -k = p[-pi=p2-p2 , 

Q = {k + k')/2 . 

The external nucleons are on shell and the following constraints hold 

= {W'' + 2^ + g')/4 , 
W-z = W-q = z-q = . 
For the Mandelstam variables, one has 

t = q , 

si = [Q^+Q-{W+z)-t/'i+m^] , 

wi = [Q''-Q-(W+z)-t/'i+m^] , 

ui = {W+z)-Q/2m , 

S2 = [Q^+Q-iW-z)-t/A+m''] , 

U2 = [Q''-Q-{W-z)-t/4.+m'] , 

V2 = {W-z)-Ql2m . 
Sometimes it is useful to write 

= (fc^ - M^)/2 + (A:'^ - M^)/2 + - t/^) , 
Q-9=(fc'^-M')/2-(fc^-M')/2. 
For free spinors, the following results hold: 

Hp) iu{p)f^ = [u(p') iu{p)f^=o, 

Hp) {W+ /) «(P)]''^ = 2m [u{p') u{p)]'-''> , 

Hp') iw- /) u{p)f^ = 2m Hp') "(p)]''^ , 

and also 

{u-yxu}'^^^ = {{W + z)x/2m[u u] - i/2m[u a^x{p' - p)" w]}''^ , 
(gVW){M m}('' = {-i/2m[u ct^a(p' -p)" u\{W + 0)^/2m}('' , 
{m7pm}'^^ = {(W^ — z)p/2m[u u] — i/2rn[u (Jup(p — p)" ^t]}'^'' , 
(gV4m''){w w}'^^ = {-i/2m[« cT.p(p' -p)" - z)'' /2m}^^'^ . 

In the CM one has 

Pi = {E;p) , p'l = iE;p') , 
P2 = {E--p) , p'2 = {E;-p'), 
W = i2E;0), 
q = (0;p' -p) , 
z = (0;p' + p) 

and the on shell condition for nucleons reads 

E'^ = m'' + qV4 + ■ 
In the CM frame, the nucleon spin functions may be expressed in terms of two component matrices 



(Al) 
(A2) 

(A3) 
(A4) 

(A5) 
(A6) 

(A7) 
(A8) 
(A9) 
(AlO) 
(All) 
(A12) 
(A13) 

(A14) 
(A15) 

(A16) 
(A17) 
(A18) 

(A19) 
(A20) 

(A21) 
(A22) 

(A23) 
(A24) 
(A25) 
(A26) 
(A27) 

(A28) 



{w(p') u{p)Y'^ = X* 

|^«(p') (^Mib'-P)" m(p)} 
= s{i) 



(i) 



/ ' \ , /' 2 . \ (P'+P) 

<Tx{p-p) + {q -^a.qxz) 



(A29) 
(A30) 



(A31) 



where s{i) = (1,-1) for i = (1,2). These results, which contain no approximations, allow one to 
write the identities 



- {[« 4"'^ [u <y,x ip' - pY - (1 - 2)} |- = 4m^ 



z 

2rn 



Asm? 

' Vlw? 
1 

"4m2 



2 2/x4\ 9 



(l+Wz^A*) 



16m4 



[m cr^xip' - pY w] [w CTi'pip' - pY w] = ^m' 



4 4 

' 16m4A4 



, , 2/-,2\ q^Z^ f2_LS , , 2/-,2\ 

+ (1 + ^ A ) .4X2 + (1 + ^ A ) 



4m4A2 ' V ' > 16m4 



where the two-component spin operators f2 were defined in section II and A^ = 4m{E + m). 



APPENDIX B: LOOP INTEGRALS 



The basic loop integrals needed in this work are 



SC 



2m^ 



[Q2+Q.(w/+2)_t/4] ' 
2m/U 



2m/L( 



[Q2+Q.(VF+^)-t/4] [Q2 + Q.(l^-^)-t/4] ' 



with 



{27rY [(Q-g/2)2-M2][(Q+g/2)2-M2] " 

All denominators are symmetric under q — > — g and results cannot contain odd powers of this 
variable. The integrals are dimensionless and have the following tensor structure: 

t{000) 



(47r)2 

i 



{nr'} , 



(47r)2 \ /i- 



{g^^g^' + g'^^g'" + g'^'g"^) nT^] } , 



(477)2 1 2m ^ ^ ^= J J ' 

/r = jiy[^ [I'l'^ nl^)] + 4^ [(w^ + ^)''(i^ + ^Y nT'^] + g"" ^fr I 



(47r)- 



{nir} 



The usual Feynman techniques for loop integration allow us to write 



with 



da {-Caf 



Jo 
8 Jo 



Po-ln(^ 



da i-Ca)' % 



-Pi+ln(^ 



p.-ln(^ 



Yf(kmn) 



2m \ 



m+n+1 n\ 



fdaafd,^-^ 
Jo Jo ^' 



2m > 






2 


2m 








2m ^ 






2 



da a I db 



-po + In ( ^ 



da a ; db b dc 





1 ,.2 



Ca = a - 1/2 



Cb = a6/2 , 

-.2 



S;:,c = -(1 - 2ab) qV4 + (1 - a6) /i' , 
Cc = abc/2 , 

The case (cs) is obtained from (sc) by making — > — The case (us) is obtained from [ss) by 
making <-> — Cc- 



APPENDIX C: OPEP ITERATION 



(Bll) 
(B12) 

(B13) 

(B14) 

(B15) 

(B16) 

(B17) 

(B18) 
(B19) 

(B20) 

(B21) 

(B22) 

(B23) 

(B24) 
(B25) 

(B26) 

(B27) 

(B28) 

(B29) 

(B30) 



The iteration of the OPEP has to be subtracted from the elastic scattering amplitude, in order 
to avoid double counting in the potential. In this work we adopt the procedure used by Partovi and 
Lomon [5] , based on a prescription developed by Blankenbecler and Sugar [41] . In this appendix we 

adapt their expressions to our rclativistic notation and also simplify some of the results. 
The iterated OPEP is contained in the box diagram, corresponding to the amplitude 



where 



"Tus — i 



4 J ^■■■J 



2mu 


(1) 


2mu 


(2) 


.u—m-' 




.s—m-' 





Evaluating this integral using the results of appendix B, one recovers the spin structure of eq.(2.6) 
with 
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n(002) , fr{000) 
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The iterated amplitude is denoted by 7^ and given by 



with 



Tit = -i 



4 

4 



— [(mm)'"'"-' (m 7i m)'^' — (m 7i m)'^^(u-u)'^-' 

r ,.2 



/c + (/b - 2 Ic) 
m 2m 



The functions are three-dimensional loop integrals, defined as 

»„3 



where = y/m? + {Q^^z/2)^ and 



)3 [(Q-q/2)2+^2][(Q + q/2)2+^2] 
The usual Feynman parametrization techniques, the representation 



(CI) 



(C2) 



(C3) 
(C4) 
(C5) 
(C6) 
(C7) 
(C8) 

(C9) 



(CIO) 



(Gil) 

(C12) 
(CIS) 
(C14) 



(C15) 



m 1 



mE de 



Eq 77 [{Q - z/2)2 + m2 + e2£2] ' 
and the tensor decomposition 

" - (47r)2 ^ " J ' 



(47r)2 [ 2rn 

(47r)2 [ " + 4m2 " + J 



(for a: = ^, B, C) yield 
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The functions 11 and n are written as 



(020) 
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where 



Pi = C\q-CbZ , 

S| = + ab{l + e^) E^ , 
T,% = T,l, + abX{l + e^)E\ 

The contribution from the OPEP cut in the functions Ilus is canceled by the integrals Ha. We 
parametrize the loop momentum in those integrals as Q = {abc W/2) = (— Cc W), and have 
[Q2 + - Pj] = Dus and write 
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The integrals lis and He can also be simplified, by adopting the new variables c and 0, defined 
by the relations e — V a'^b'^c'^ — ah cos 6 /Vab, Q = E\/a'^h'^c^ — ah sinO. Performing the angular 
integrations, we have 
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The results presented so far in this appendix correspond just to a reorganization of those obtained 
by Partovi and Lomori [5] . They may be further simplified by noting that 

lB=i [{ 
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The integrals /(•••) can be performed analytically and we have 
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where Ha is the function given in eq.(7.5). 
Thus, 
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The results presented in eqs.(C3)-(C8), (C20)-(C23), (C33)-(C35), and (C43)-(C47) allow one 
to write 
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where the integrals Xl^eg = ^it — ^us are regular and given by 
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APPENDIX D: FULL RESULTS 

In this appendix we list the results for the amplitudes that enter eq.(2.15), obtained by reading the 
diagrams of fig.6 and representing loop integrals by means of the functions displayed in appendices 
B and C. 

family 1 (diagrams a+b+c+d+e+f) 
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APPENDIX E: RELATIONS AMONG INTEGRALS 



We display here the relations among integrals needed for the chiral expansion of the potential. 

The derivation of these relations is based on the fact that the numerators of some integrands can 
be simplified. For instance, a result for /^^ may be obtained through 
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where the ellipsis indicates that short range contributions were discarded. The combination of both 
results produces 
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Other two relations involving 113°°°^ and U.^°g^ are obtained by deriving eq.(B20) and eq.(C55) 
with respect to n, 
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APPENDIX F: INTERMEDIATE RESULTS 

The results presented here for the TPEP were obtained by using the relations among integrals of 
the preceding appendix into the full expressions of appendix D. In this procedure we just neglected 
short range integrals and both sets of equations are equivalent for distances larger than 1 fm. In 
family 3, we did not keep contributions larger than 0{q*), in order to avoid unnecessarily long 
equations. 

family 1 (diagrams a+b+c+d+e+f) 
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family 2 (diagrams g+h+i+j) 

• = - fi - '^Z'^') - - 2-] ' , 

and 



+ M ^ - ) ((l-t/4M^) + 2 - i/4,x^ 



4/^2 



3m I /I 



(47r)4 



~ (4^)4 



^((i-VV)n(r + ^(i-V2M^)n 



(001) 



(g)_ _ _/AnV4/^ 



^ ((i-t/4M^) nr^ + ^ (i-t/2M^) nr^' - ^) 



2 2/0-^2 2 

(w)- ^ _l£jri_/jfU_ 9_ tt{002) 
(47r)4 m2 



(000) _ JL (i-t/2M') nl°°i' 



^(z)- _ fJ- ■in n-n: y „(002) 



2 2 I £2 2 



y 3m y m^ ' 



family 3 (diagrams k+l+m+n+o 

,3/o^J-2 „2 



(000) 



+1 (4o + Vm' 4i) <5i + i (5^)' n; 



(000) 
cc ) 



(47r)2 m2 

)+ _ mV2/^ 



{ [4o + (V/x^) St.] nT'' + 1 sto (i-t/V) nr^ } 



'-DB 



(47r)2 rri 



^ { [5o+o + (Vm7 4] ni""^) - 5to (1- V V) nr ^ } , 



~ (47r)2 ^ 3 U j iicc 



and 



r('")- 



^os --74^/^o-o|3(^-^)(l-t/V)n, 



2n T-r(OOO) 



■^DB ^i3S ' 



(i-t/4M^)nir + ^(i-V2M^)n(r' 



APPENDIX G: RELATIVISTIC EXPANSIONS 

In section VII we have discussed the relativistic expansion of the function ^{t) derived by Becher 
and Leutwyler, which does not coincide with the usual heavy baryon expansion. In this appendix 
wc show how their results can be used to produce relativistic expansions for box and crossed box 
integrals. 

The triangle, crossed box, and regularized box integrals given, respectively, by eqs.(B17), (B19) 
and (C54) can be written as 



nir = - ['dcuTHMs.), 

Jo 

/oo 



where nf°"\M) is a generalized triangle integral, given by 

n + l ^1 „i 



and the denominator D{A4) is 
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D{M) = M'^a'^b'^ - a{l~a){l~b) q'^ + (1 - ab) . 

When M. = m, one recovers the triangle integral defined in eq.(B17). On the other hand, the values 
M'^ = {W^ + q^ + c^ z'^)/A and = (c^ + + z'^)/A yield eqs.(Gl) and (G2). 

Performing explicitly the 6 integration in eq.(G4), we obtain the generalization of eq.(E2), which 
reads 
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with 



JL 



In 



AT 



2 +2 



M 



4M^ 



tan 



1 / M ll-n/2M\ 



M V 1 + M/2A1 y 



In all cases is a large parameter and we can use the relativistic expansion of 7(t), which is 
related to our triangle integral by Ilic""'' = — 2m/Li(47r)^ 7(^)- We have 
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with Da and 11^ given by eqs.(7.5) and (7.6). Recalling that H, 
into eqs.(Gl) and (G2), we obtain 
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and inserting these results 
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In the chiral limit /^ — » 0, we have 
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and, using cqs.(E7), (E8), (E9), and (E17), and (E12), (E13), (E14), and (E18), one finds the 
following relationships valid in that limit 
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These results may also be combined with those presented in appendix E, in order to produce 
relativistic 0{q^) expansions for box and crossed box integrals. Eqs.(E8), (G13), and (E2) yield 



n 



(000) 



- - 1 + ^— + 



and, using eqs.(E17) and (E7), one has 



n 
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Recalling that nis""^ = IIx and using the results of section VII, we find the heavy baryon 
expansion 

,2 



m (l-t/V) 4m2 

where the ellipsis represents polynomials in t. 

For the box integrals we evaluate eq.(G8) directly and obtain 
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Comparing with eq.(E14) and using eq.(E2), we find 
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where lib = li^e^-' . Finally, evaluating e(j.(E18), we have 
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